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By analysizing the global properties of the equations for 
general relativistic collapse which reflect the rate of change 
of the circumference radius R with the proper radius I, we 
show that (i) no "trapped surface" is ever formed and (ii) the 
gravitational mass of the fluid M(r,t) — > as R —* 0. This 
result is absolutely independent on the equation of state of 
the fluid or any other details. Thus, for continued collapse, 
the total energy output measured by a distant observer, is 
Q — > A/jC 2 , the original gravitational mass. 

PACS numbers: 04.20.Dw, 04.70. Bw 



The problem of gravitational collapse has been de- 
scribed as "the greatest crisis in physics of all time" 
(pp.1196). Though the theoretical framework for 
studying this problem is fairly well developed (2^, in 
practice it is impossible to make much headway without 
making a number of simplifying assumptions because of 
our inability to self consistently handle: (i) the equation 
of state (EOS) of matter at arbitrarily high density and 
temperature, (ii) the opacity of nuclear matter at such 
likely unknown extreme conditions, (iii) the associated 
radiation transfer problem and all other highly nonlinear 
and coupled partial differential equations. As far as ex- 
act solution is concerned, it is possible only if the collaps- 
ing fluid is considered to be a dust with pressure p = 
everywhere including the center even for arbitrary high 
density. Further, it should be considered homogeneous, 
and the pioneering solution to this effect was obtained by 
Oppenheimer and Snyder 0](OS) which suggested that 
the spherical dust ball collapses to a point to a state of 
infinite proper density, i.e., R — * 0, p — > oo in a proper 
time (by fiat, pressure still remains zero): 
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where p(0) is the initial mass-energy density of the fluid 
when it is assumed to be at rest. A distant inertial ob- 
server 5*00, however can not witness this approach to a 
physical singularity because as R — > Rg — 2GM i jc 2 , 
where M, is the initial gravitational mass of the dust, 
the photons emitted from the suface of the collapsing 
fluid are infinitely redshifted indicating the formation of 
an "event horizon" and a subsequent black hole (BH). 
Attempts for numerical studies fa] to study the collapse 
of a physical fluid can progress only by presuming that 



the energy emitted during the process Q <C Af^c 2 . The 
said assumption tantamounts to the condition that the 
eventual temperature of the nuclear matter T <C lGeV, 
the typical nucleon mass. When this is so, one can use a 
T = degenerate EOS or can incorporate a finite T EOS 
which incorporates a first order correction term cx T 2 . A 
general overall self-consistent numerical treatment of the 
problem which does not impose explicit or tacit smallness 
conditions on Q and T is implausible. For a dust, ther- 
modynamics would tell that a p — condition implies 
that the internal energy density e = and also T = 0. 
Consequently, the instantaneous value of Q = 0, which in 
turn implies the instantaneous gravitational mass of the 
body M = Mi. The presumption that Q <C Af^c 2 leads to 
a more or less similar condition, that M « Mi, and natu- 
rally numerical computations reveal the formation of an 
event horizon enclosing a mass M g = Mf « Mi within a 
finite proper time. However, when the condition that for 
all times Q <C Af^c 2 is not imposed the running value of 
M = Mi — Q jc? < Mi , consequently, the would be value 
of R g will be smaller than the initially suspected value of 
i? g i leading to value of Q higher than the previous value. 
This may in turn lead to a value of M lower than its 
previously expected value. And, in general, it can not be 
predicted whether this apparently runaway process con- 
verges to a certain finite value of Mg in a finite time or 
not, where M g is the mass at R = R g . In this situation, 
we would like to explore the limiting properties of the 
equations for collapse for a perfect fluid without making 
any conceptual assumption. The matter part of the fluid 
has the energy momentum tensor 



= (p + p)u^u v +pg< 
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where u M is the fluid 4- velocity, and p = po + e, with 
Po = mn, m is the rest mass of a baryon as measured by 
Sqo and n is the proper number density of the baryons. 
Now we have taken c = 1 . The radiation part of the fluid 
is described by 



= qk^k v 



(3) 



where q is both the energy density and the radiation flux 
in the proper frame and k^ = (1; 1, 0, 0) is a null geodesic 
vector H so that fc M fc M = 0. Irrespective of whether a 
"trapped surface" is formed or not, the interior of the 
collapsing body may be described by the Schwarzschild 
metric §-§§: 
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R(r,t) 2 (dtf 2 + sin 2 tfdip 2 ) (4) 



Here r is a comoving coordinate, so that the total num- 
ber of baryons N(r) enclosed within the sphere of con- 
stant coordinate radius r is constant and independent of t 
and R{r, t) is the circumference coordinate. The explicit 
forms of e^ and e A are not known, but they are nonsin- 
gular for any finite R, and computational techinques are 
supposed to self-consistently evaluate them. The exterior 
region must be described in terms of Vaidya's retarded 
time u , which turns out to be the proper time of Soo 



ds 2 = -e^du 2 



2eVdu dr + R 2 {d$ 2 + sin 2 tfdtp 2 ) 



(5) 



where e 2 ^ = (1 — 2GM(u)/R). The comoving time and 
radial derivative operations are respectively defined by 



A = e~+ 



d 



A = e 



-A 



dl 



(6) 



Here dr — e^dt is an element of proper time at constant 
r, and dl — e x dr is an element of proper radial distance 
between two baryon-shells r and r+dr, ■& and tp are frozen 
in spherical symmetry. The circumference velocity of the 
fluid measured at constant r in terms of proper time 



dR 



DtR = U = 

dr 



(7) 



Here D t and D r are to be treated as equivalents of 
"convective derivative" following the trajectory of the 
fluid parcel used in classical hydrodynamics. The cen- 
tral global character of the collapse process, valid in both 
metrics, is manifest through 

p ]j ft d(circumference)/2n dR 
d(proper distance) dl 

And the most important global constraint on the field 
equations is that the relationship 



T 2 = 1 - 



2GM 



R 



U 2 
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is preserved for all times ]2J,|3|]. The actual equations for 
containing the details of the collapse process, and valid 
in metric (4), are [|1 



D T M = AttR 2 [T(p + q) + Uq] 



D t M = -4irR 2 pU - L(U + T) 



D t U 



dp 

p + p \dR 



M + ATrR 3 (p + q) 
R 2 



(10) 



(11) 



(12) 



AT 



UT f dp 



p + p \dR 



L 
R 



where the comoving luminosity is 
L = AirR 2 q 



(13) 



(14) 



For motionless matter in flat space (S'oo), T = 1. In a 
curved space V < 1 because of the gravitational field. 
And a value of T < would signify that the spacetime is 
pinched off from the rest of the universe ||, and this is 
what is believed to happen in the final stages of gravi- 
tational collapse. However nobody has ever hinted what 
could be order of magnitude of the (-ve) value of final Tf 
even for the idealized case of a dust collapse. The popu- 
lar conviction that even a physical fluid must eventually 
behave like a O-S dust is based on certain properties of 
the radiationless collapse equations j^]: 



D r M = AttR 2 pT 



D t M = -AirR 2 pU 



D t U = N£ - 
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dp 
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The traditional argument favouring BH formation goes 
like this ||: For collapse, U < 0, and from eq. (16) it 
follows that M(r,t) increases monotonically within the 
fluid; of course, it remains constant, for adiabatic col- 
lapse, for R > Rq, the outer boundary. Because of this, 
the M/R 2 term in cq. (17) increases more rapidly than 
than it would in the Newtonian theory where it remains 
constant. Then eq.(18) would suggest that T decreases 
monotonically instead of remaining constant in a Newto- 
nian case. Both of these relativistic effects are supposed 
to cause the collapse process monotonic and accelerate 
faster than the Newtonian case S. Here note that un- 
less we presume that dp/dR can be positive, a negative 
value of r in eq.(18) would demand that T increases with 
time implying T would tend to be back to zero. In any 
case when we restore q, or in other words, when we do 
away with the assumption Q <C Mi, this entire line of 
argument collapses and nothing can be predicted about 
the nature of evolution of the collapse process. In this 
situation, we discuss below the global properties of these 
equations. Following the discussion by Landau & Lif- 
shitz J! (pp.250), the ordinary radial velocity of the fluid 
as measured in a local rest frame is 



A dr 
*~di 



dl_ 



(19) 
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By principle of equivalence (POE), v < c. By com- 
bining eq. (7), (8), and (20), we find 



U = Tv 



(20) 



The above relationship does not depend on the signa- 
ture or form of the metric coefficients and, as mentioned 
earlier, all equations involving interior Sch. metric coef- 
ficients remain vaild even if a trapped surface is formed 
||. Further, eq. (20) must not be confused with u r = yv, 
where u r is the radial component of the 4-velocity mea- 
sured in the local inertial frame (in the comoving La- 
grangian frame, u r = 0) and 7 = (1 — v 2 )^ 1 / 2 is the 
Lorentz Factor of the fluid in the same inertial frame. 
By combining eq. (9) and (19), we obtain 



1 



2GM 
R 



(21) 



Since the LHS of the foregoing eq. is positive (even if 
r were negative), so will be the RHS, and therefore, by 
simply using POE, we find that 



2GM 



R 



< 1 



(22) 



Recall that, for spherical collapse the condition for for- 
mation of a trapped surface, a surface of no-return is 
2GM/R > 1 gQ. Thus, in a most general manner, 
we find that trapped surfaces are not allowed by GTR 
contrary to the assumtion made in a large number of pre- 
vious works staring from the sixties. It is the assumption 
(1) about the occurrence of the trapped surfaces, along- 
with, (2) the causality conditions requiring the nonexis- 
tence of closed timelike lines, (3) the strong energy con- 
dition (T M „ — 0.5g M „)u M 'u ly > 0, and (4) a certain gener- 
ality condition on the Riemann-Christoffel tensor, that 
laid the foundation of the famous singularity theorems 
of Hawking, Penrose, Geroch and others |Io[| . So now 
they are invalidated because trapped surfaces are not al- 
lowed by GTR. We quickly recall the relevant comments 
by Misner et al. Q (pp.935): "All the conditions except 
the trapped surface seem eminently reasonable for any 
physically realistic space time" . Thus, the mathematical 
and moral foundation for the inevitability of occurrence 
of singularities in GTR are eliminated. However, this 
does not necessarily mean that singularities must not oc- 
cur, and, we will probe this question in an accompanying 
paper | p0[ . Any EOS will support only a certain finite 
maximum Mi and therefore given sufficiently large Mi, 
continued collapse is a genuine possibility; and if indeed 
R — > 0, then in order to satisfy the foregoing ineqality, 
we must have 



M f -> 0; 



R 



f 







(23) 



Remember here that the quantity Mq = mN (which is 
the baryonic mass of the star, if there are no antibaryons) 



is conserved as Mf — > 0. How do recoincile our general 
result Mf = 0, with the result of O-S collapse which 
suggests that as the dust ball with an initial condition 
R = Ro, M(r o ,0) — Mi collapses from a state of rest, 
i.e., v = v(r, 0) = and U — U(r,0) = 0, it forms a 
BH with Mf = Mi in a finite proper time? To find the 
p(0) occuring in eq. (1) in a truly self-consistent way, 
we must fix it from the collapse equations themselves, 
rather than artificially and by hand. When we put this 
U = q = condition in eq.(12), we obtain the so-called 
Oppenheimer-Volkov (OV) equation pT[. 



dp P + p(0) 

dR ~ ~ R(R — 2GM) 



(AirpR 3 + 2GM) 



(24) 



By feeding the p = EOS in the OV equation it may 
be verified that whether the fluid is considered homo- 
geneous or not the solutions could be of two types: 
(i)p(0) = 0, R = 00, Mi = finite and (ii) p(0) = 
0, Rq = finite, Mi = 0. In the former case, the fluid 
boundary requires infinite r for attaining a finite R or 
R g , implying that the collapse process hardly ever takes 
off: 



r(R) = 



Ro 

c 



Ro 

R e 



1/2 



R 



1/2' 



(25) 



and there is no question of formation of a BH. And in 
the latter case, our proof that Mf = is vindicated be- 
cause Mf = Mi = 0. Even in this case, eq. (1) yields 
r = 00! These boundary conditions are fudged in 
all theoretical and numerical discussions by pretending 
that a self-supporting fluid with finite values of p and 
Rq, and p / 0, everywhere, may approximately be con- 
sidered to suddenly attain the dust-like dynamics. Even 
if such a dust-like collapse took place, it can be verified 
from eq.(18) (by evolving from U — condition) that, 
far from assuming any negative value, the value of T re- 
mains fixed to the initial value Tq, which is unity, and 
which is the appropriate value of T for a motionless (no 
kinetic energy) and matter free (thus gravityfree) region. 
It might be thought that if the numerical calculations 
study the collapse of a supermassive star, the horizon 
may be formed at very small densities ~ 10 _4 g/cm 3 , 
and the entire collapse physics including the EOS can be 
handled accurately. This idea would not be correct be- 
cause, the expectation that horizon would be formed at 
R g « 2GM{ is based on the presumption that Q <C Mj. 
And in any case our central result manifest by eq. (22) 
is an exact one. In the context of a numerical neutron 
star collapse, we will discuss about an ansatz for handling 
the thremodynamics of the problem which may reveal a 
value of Q much higher than previously reported values, 
although the eventual computation must break down as 
T CO m keeps on increasing without any upper limit p3[ . 

It is widely believed that Chandrasekhar's discovery 
that White Dwarfs (WD) can have a maximum mass set 
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the stage for having a gravitational singular state with 
finite mass. The hydrostatic equilibrium of WDs can be 
approximately described by Newtonian polytropes Jl3] 
for which one has R oc p^ n ^/ 2 " ) where p c is the cen- 
tral density of the polytrope having an index n. It shows 
that, for a singular state i.e., for p c — > oo, one must have 
R — > for n > 1; and Chandrasekhar's limiting WD 
indeed has a zero radius. On the other hand, the mass 
of the configuration M oc p]? n ^ 2n . And unless n = 3, 
M — > oo for the singular state. However, for ultrarel- 
ativistic electrons, the EOS is p — > e/3 and the corre- 
sponding n — > 3 and one obtains a finite value of - 
the Chandrasekhar mass. Now when we apply the theory 
of polytropes for a case where the pressure is supplied by 
the baryons and not only by electrons, we must consider 
GTR polytropes p5[ . It can be easily verified from eq. 
(2.24) of ref.(15) that in the limit p c — > oo, the scale size 
of GTR polytropes A^ 1 -> 0. Further eq.(2.15) and (4.7) 
of ref. (15) tell that M oc K n ' 2 oc p~ 1/2 -» for p c -> oo. 
Thus, a proper GTR extension of Chandrasekhar's work 
would not lead to a BH of finite mass, but, on the other 
hand, to a singular state with M — > 0. Physically, the 
M = state may result when the negative gravitational 
energy exactly cancels the internal energy, the baryonic 
mass energy Mq and any other energy, and which is pos- 
sible in the limit p — ► oo and p — > oo. In fact, Zeldovich 
and Novikov jl2| discussed the possibility of having an ul- 
tracompact configuration of degenerate fermions obeying 
the EOS p — e/3 with M — > 0; we have verified that their 
configuration actually corresponds to R — > R g — > 0. In 
fact, they mentioned the possibility of having a machine 
for which Q — ► M^c 2 . In a different context, it has been 
argued that naked singularities produced in spherical col- 
lapse must have Mf = |L6]]. Also, it was discussed long 
ago that spherical gravitational collapse should come 
to a decisive end with = M* = 0, and, in fact, this 
understanding was formulated as a "Theorem" |l7j ]: 

"THEOREM 23. Provided that matter does not un- 
dergo collapse at the microscopic level at any stage of 
compression, then, -regardless of all features of the equa- 
tion of state - there exists for each fixed number of 
baryons A a "gravitationally collapsed configuration" , in 
which the mass-energy M* as sensed externally is zero." 
We have proved this and more for spherical symmetry, 
our result is absolutely general and has got nothing to 
with the exact EOS and radiation transport equations. 
On the other hand, it depends only on the global prop- 
erties of the GTR collapse equations and the Principle of 
Equivalence. We have learnt that GRB990123 released 
an energy 3 x 10 54 ergs of energy in 7-rays alone under 
condition of isotropy. Given the fact that the relativis- 
tic blast wave radiating the gamma-rays may not have 
an efficiency greater than 10-20% in converting the blast 
wave energy into gamma rays, the total electromagnetic 
energy released in the process could easily be Qfb ~ 10 55 



erg. Further, the energy output in us could be or rather 
should be > Qfb , the net energy released in the process 
could be Q ~ few 10 55 erg. In case, there is weak beam- 
ing, the total energy released may still be Q < 10 55 erg. 
physical mechanism, other than what has been propsed 
here, can satisfactorily explain the origin of such unprece- 
dented energy. Such energy release is hardly possible if 
trapped surfaces would form. 
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